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ABSTRACT 

We present new results on the star formation history of the dwarf irregular galaxy NGC 1569. The 
data were obtained with Hubble Space Telescope's NICMOS/NIC2 in the F110W (J) and F160W (H) 
near-infrared (NIR) filters and interpreted with the synthetic color-magnitude diagram method. The 
galaxy experienced a complex star formation (SF) activity. The best fit to the data is found assuming 
three episodes of activity in the last 1 — 2 Gyr. The most recent and strong episode constrained by 
these NIR data started ~ 3.7 x 10 7 yr ago and ended ~ 1.3 x 10 7 yr ago, although we cannot exclude 
that up to three SF episodes occurred in this time interval. The average star-formation rate (SFR) of 
the episode is ~ 3.2 M yr -1 kpc -2 , in agreement with literature data. A previous episode produced 
stars between ~ 1.5 x 10 8 yr and ~4x 10 7 yr ago, with a mean SFR about 2/3 lower than the mean 
SFR of the youngest episode. An older SF episode occurred about 1 x 10 9 yr ago. All these SFRs 
are 2 — 3 orders of magnitude higher than those derived for late-type dwarfs of the Local Group. In 
all cases an initial mass function similar to Salpeter's allows for a good reproduction of the data, but 
we cannot exclude flatter mass functions. These results have been obtained adopting a distance of 
2.2 Mpc and a reddening E(B-V)=0.56. A larger distance would require younger episodes and higher 
SFRs. We have explored some possible scenarios using the astrated mass in the best fit model, in 
order to constrain the past star formation history. We cannot rule out a low rate in the past SF but 
we can safely conclude that the last 1 — 2 Gyr have been peculiar. 

Subject headings: galaxies: evolution — galaxies: individual: NGC 1569 — galaxies: irregular - 
galaxies: dwarf — galaxies: stellar populations 



1. INTRODUCTION 

In the last decade the field of galaxy evolution has 
gained substantial momentum, owing to the wealth of 
high-quality multi-wavelength surveys from the main 
ground and space-based telescopes and sophisticated the- 
oretical models. Nonetheless, the current knowledge of 
the major physical mechanisms for galaxy formation and 
evolution is still affected by serious uncertainties. Tra- 
ditionally, star formation (SF) in cosmology and galaxy 
formation studies is naively described by simple univer- 
sal laws, and the stellar initial mass function (IMF) is 
assumed to be constant in time and space. The only way 
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to obtain detailed information on these crucial evolution 
parameters is to study the stellar populations of nearby 
galaxies and derive the SF history (SFH) over cosmologi- 
cal timescales. One of the fundamental tools is the color- 
magnitude diagram (CMD) of the galaxies resolved into 
single stars. Since all the CMD features are related to the 
evolutionary state of the stellar populations, we can un- 
derstand the SFH of a certain region within a galaxy by 
interpreting its CMD in terms of stellar evolution theory 
with the method of the synthetic CMDs. In theory, from 
the CMD we can infer how many bursts have occurred, 
their duration and star formation ra t e (SFR), and the ini- 
tial mass function ijTosi et alJll991l: IGreggio et alJll99S|) 
for any star-formation law and metallicity. In practice, 
observational errors and theoretical uncertainties impose 
severe restrictions. The requirement of resolving a galaxy 
into stars limits the exploration of the SFH to nearby 
galaxies. At present, even the largest telescopes only 
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permit such studies for galaxies with distance < 15 — 
20 Mpc, i.e., in the local universe. In this limited space 
volume we find many spirals and a plethora of both late 
and early dwarf galaxies. 

Dwarf irregular galaxies (dlrrs) and blue compact 
galaxies ( BCDs) play a major r ole in the field of galaxy 
evolution ijGrebell Il997t iMated 119911 . With their high 
gas content and lo w metallicity, they could be similar to 
primeval galaxies ijlzotov et al.lll997|) and the best possi- 
ble sites for determining the primordial He a bundance 
Ijlzotov fc Thuanl 120041 lOlive fc Skillmanl 12004ft . It has 
been suggested that the excess of number counts of faint 
blue galaxies at redshift z ~ 1 may be ascri bed to dwar f 
galaxies undergoing their first burst of SF lpIhllT997|) . 
Dwarf galaxies are important in galaxy formation mod- 
els since in hierarchical models they are assumed to be 
the building blocks from which other types of galaxies 
form through merging. 

Key questions on the evolution of late-type dwarfs are: 
did all the dlrrs and BCDs experience SF in the past? If 
so, what are the properties of the SF episodes? Several 
stellar population studies have suggested that dlrrs are 
characterized by a gasping (episodes of moderate activity 
separated by shor t quiescent phases) SF with no evidence 
of extended gaps jTqsi et al.|199lHSchulte -Ladbeck et al.l 
120011 iDolnhin et all 12003 iSkillman et al] 120031) . Thek 
SFR is not strong enough to explain the faint galaxy 
counts excess. 

NGC 1569 is a dlrr near the Local Group. Its 
intrinsic dista nce modulus has been estimated to be 
26.71 ± 0.60 (|Arp fc Sandagdll985l Iisrae]|ll988l: iWallerl 
1 19911). corresponding to a dist a nce of 2.2 ± 0.6 Mpc. 
lo'Connell. Gallagher fc Hunterl l|1994D obtained a dis- 
tance modulus o f 27.0 ± 0.5 corresponding to 2.5 
± 0.5 Mpc, while iKa rachent sev. Tikhonov fc Sazonoval 
lfT99l and IRicher fc McCalll (|1995l) found 1 .8 + 

0.4 Mpc and ~ 1.7 Mp c, respectively. Recently, 

iMakarova fe Karachentsevll)2003|) derived a new distance 
to NGC 1569 based on the I magnitude of the tip of the 
red giant branch (RGB). Their data are consistent with 
two possible distances: 1.95 ± 0.2 or 2.8 ± 0.2 Mpc. 
Herea fter we adopt a distance of 2.2 ± 0.6 Mpc ({Israeli 
1988), which encompasses the wide range of distances 
proposed in the literature. 

The low Galactic latitude (b = 11 2°) of NGC 1569 
implie s significant Galactic reddening. iBurstein fc Hcilcs 
l|1984j) estimated a value of E(B - V) = 0.51 for th e fore- 
groun d extinction using the HI column density. Ilsraell 
l|1988() found a total extinction of E(B - V) = 0.56±0.10 
from the UV color-color diagram. Furthermore, many 
authors have suggested t he existence of a strong, d i fferen- 
tial, internal reddening. iGonzalez-Delgado et al.l l)1997f) 
found a total reddening of E(B - V) = 0.67 ± 0.02 in 
the direction of the super star clusters (SSCs), where the 
component of the intrinsic red dening is E(B — V) = 0.11. 
iDevost. Rov fc Drissenl l)1997j) estimated the Galactic ex- 
tinction to be E(B — V) = 0.50 and the intrinsic mean ex- 
tinction E(B — V) = 0.20, using optical emission spectra 
of the ionized gas in NGC 1569. iKobulnickv fc Skillmanl 
l)1997[) found a variation of the total reddening in the 
galaxy from E(B - V) ~ 0.51 to E(B - V) ~ 0.77 with 
an average value E(B - V) ~ 0.63 (with R = 3.1). They 
found it reasonable to assume a Galactic foreground of 
E(B — V) = 0.51 and ascribe the difference to the in- 



ternal reddening. In IKobulnickv fc Skillmanl l)1997fl . the 
extinction for the area discussed in the present paper 
is E(B -V) ~ 0.58, except for the SSCs and HII re- 
gions. lOriglia et alJ l)2001|) derived a foreground ex- 
tinction of E(B — V) = 0.55 ± 0.05 and an intrinsic 
reddening of E(B — V) = 0.15 ± 0.05 from UV spec- 
tra. The simulations p resented in this paper assume 
E(B - V) = 0.56 ± 0.10 ll!sraeillT98^ . 

The mean oxygen abundance of NGC 1569 found in 
the literature is 12 + log(0/H) ~ 8.2 ± 0.2 dex. With 
an assumed [0/Fe]= 0.0, the average met allicity is then 
0.25 Zq corresponding to Z ~ 0.004 (see iGreggio et all 
1998, and references therein) if Z© = 0.02 2 . 

This rather exceptional dwarf galaxy is known to con- 
tain two SSCs, SSC-A and SS C-B, with SSC-A hav- 
ing two comp onents, Al and A2 l)De Marchi et al.lll997t 
lOriglia et al1l200l|) . Many HII regions a re observed i n 
those zones where SF is currently active (|Wallerlll991[) . 
Nevertheless, NGC 1569 is considered a "post" star- 
forming galaxy with a total dynamical mass of M^yn ~ 
3.3 x 10 8 Mq (|IsraelHl988j) . one-third of which is neu- 
tral h ydrogen gas (M H i ~ 1-3 xlO 8 M^. IStil fc IsTadl 
2002). This galaxy is a gas-rich system in a relatively 
early stage of chemical evolution. Many filaments are 
found at different wavelengths, and there is a strong spa- 
tial correlation bet ween the extended X -ray emission and 
the H Q filaments (|Martin et al.ll2002fl . The HI emis- 
sion of NGC 1569 shows a dense, clumpy ridge distri- 
bution surrounded by a more extended diffuse neutral 
hydroge n. There are discre te features su ch as arms and 
bridge s (jStil fc Israeli l2002|) . In the IR, iLisenfeld et"ai] 
( 2002) found both large an d very small grains e xposed 
to a strong radiation field. Gall iano et alJ l)2003l) found 
different dust properties in NGC 1569 (most of the grains 
are small) compared to other more metal rich galaxies. 
These results are consistent with the presence of shocks 
having dramatic effects on the dust. In the optical, deep 
H a images show filaments and arcs of warm ionized gas. 
These filaments have a length of several kiloparsecs with 
typical velocities around 50 — 100 km s _1 . The high- 
est velocities, up to 200 km s _1 , have been attributed 
to ex panding sup erbubbles with dynamical ages of 10 
Myr (jWallerlll991|) generated by a strong SF episode. At 
high energy, X-ray observations reveal that ~ 60% of the 
soft X-ray emission in NGC 1569 is centered between 
the two SSCs. Its origin is attributed to thermal emis- 
sion from the h ot gas of a galactic sup erwind emanating 
from the disk l)Heckman et alJll995l) . Chandra obser- 
vations find many high-mass X-ray binary systems and 
large inhomogeneities in the metal abundances of the in- 
terstellar medium (I SM) with ranges from 0.1Z & to IZq 
(jMartin et al.ll2002j) . This is evidence that the ISM has 
been affected by n umerous supernova e xplosions from a 
recent SF episode llMartin et al.M2002D. Hydro dynami- 
cal simulations of D' Ercole fc Brighentil l)1999t) demon- 
strated that a galactic wind triggered only by SNII ex- 
plosions (i.e. without stripping or other environmental 
effects) would not be very effective in removing the ISM 
in NGC 1569. In fact, even if most of the metal-rich stel- 
lar ejecta can be efficiently expelled from the galaxy and 

2 In the last years, the solar abundances have undergone major 
revisions (Asplund et al. 2004). The current estimate of the solar 
metallicity is Z0 inew = 0.013 
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dispersed in the intergalactic medium, cold and dense gas 
replenishes the central region of the galaxy a few hundred 
million years after the starburst. 

The SFH of NGC 1569 has most recently been studied 
by two groups with Hubble Space Telescope (HST ) data . 
In optical bands (V and I) , iVallenari fc Bomansl l|1996f) 
found a global SF episode from 100 Myr ago to 4 Myr 
ago and hints of an older episode of SF from 1.5 Gyr ago 
to 150 Myr ago. They ruled out the existence of long 
quies cent phases in the la st 1.5 Gyr. Based on B and V 
data, iGreggio et al.l l)1998l hereafter G98) found a global 
SF episode of > 0.1 Gyr duration, ending ^5 — 10 Myr 
ago. During the burst, the SFR was approximately con- 
stant, and, if quiescent periods occurred, they lasted less 
then ~ 10 Myr. The derived SFR is very high and equal 
to 3, 1 or 0.5 M yr _1 for a single slope initial mass func- 
tion with an exponent a = 3, 2.6 or 2.35, respectively, 
and a mass r ange from 0.1 to 12 M Q . From HST NIC- 
MOS images. lAloisi et alJ (|20Qll hereafter AL01) found 
that young stars (M > 8 M , age < 50 Myr) are more 
clustered around the two SSCs, while the intermediate- 
age objects (1.9 Mq < M < 8 M Q , 50 Myr < age < 1 
Gyr) have an almost uniform distribution, and the old 
stars (M < 1.9 M Q , age > 1 Gyr) are located at the 
outskirts of the starbursting region. 

All these studies show that NGC 1569 experienced a 
strong SFR in recent epochs. With a total mass of ~ 
3.3 x 10 8 Mq, the system cannot have sustained such a 
strong activity over a Hubble time, since it would have 
consumed all the gas in less than 1 Gyr (G98). The 
question then is whether the SF in NGC 1569 has started 
only in relatively recent epochs or has lasted for several 
gigayears but at much lower rates. Our goal is to obtain 
a better understanding of the past history of this galaxy. 
The best way to study old stellar populations, at least in 
principle, is to observe in the IR bands where the RGB 
phase is more evident, and where reddening effects are 
not as severe as in optical bands. 

In this work , we derive the SFH with the synthetic 
CMD method |Tosi et alJ l|199lD and G98] using the data 
of AL01 in the NIR filters. In Section 2 we summarize 
the data and the observed CMD. In Section 3 we briefly 
describe the CMD method and the artificial star exper- 
iment in order to derive completeness and photometric 
errors. The results are presented in Section 4. They are 
discussed and compared with the literature in Section 5. 

2. FROM THE OBSERVATIONS TO THE CMD 

We briefly summarize the main characteristics of the 
observations and photometric reductions; for further de- 
tails, see AL01. The observations of NGC 1569 were 
performed in 1998 with the NICMOS/NIC2 camera on 
board of HST in the F110W and F160W filters. The im- 
ages cover the crowded central region of the galaxy, where 
SSC-A, SSC-B, cluster #30, and part of an extended HII 
region complex are present. The NIC2 camera has a field 
of view of 19.2 x 19.2 , corresponding to 205 x 205 pc 2 
at a distance of 2.2 Mpc. The observations were taken 
with dithering. The resulting "drizzled" images have a 
point-spread function (PSF) FWHM of 3.0 and 3.9 pixels 
in F110W and F160W, respectively, and a pixel size of 

3 Throughout this paper we parameterize the IMF as <p(m) ex 
m~ a . In this notation, Salpeter's slope is a = 2.35 
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Fig. 1. — Top panel: observed mpnow vs. mpnow — mFi60VK 
color-magnitude diagram. Middle panel: Padova tracks of metal- 
licity Z = 0.004 for 1, 1.5, 1.7, 2, 2.5, 3, 4, 5, 6, 9, 15, 20, 30 M 
are shown as dashed lines. The boxes described in the text are 
plotted with solid lines: BOX1 (18.5 - 20.5 and 1.0 - 1.4) and 
BOX2 (20.5 - 22.5 and 0.9 - 1.3). Bottom panel: the luminosity 
function with a step of 0.25 mag. 



0.0375 , which corresponds to ~ 0.4 pc. 

The photometric reduction has been performed by 
AL01 using DAOPHOT via PSF fitting. PSF-fitting 
photometry is not the ideal procedure for drizzled im- 
ages, because drizzle creates a position-dependent PSF 
smeari ng and produces correlat ed noise between adjacent 
pixels l|Fruchter fc Hookll2002|) . However PSF-fitting is 
the only viable technique in a crowded field like the cen- 
tral region of NGC 1569. AL01 used the deepest image 
(F110W) to locate the stars and used their coordinates 
to force photometry in the F160W image. This implies 
that all of the objects detected in the F160W image are 
also present in F110W but the opposite is not necessar- 
ily true. Then, they selected the objects by adopting the 
DAOPHOT parameters x 2 < 1-5 an d sharpness < 0.4 
and obtained the final catalog of 3177 stars. The objects 
rejected by the selection criteria described above are typ- 
ically the brightest objects in the SSCs, smaller star clus- 
ters, extended objects, and blends. The star magnitudes 
were calibrated in the HST VEGAMAG system. 

The CMD and luminosity function (LF) for NGC 1569 
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Fig. 2. — The results of the artificial star experiments. The left 
(right) panel shows the ASE results for the F110W (F160W) filter. 
The completeness functions are in the top panels. The solid lines 
refer to the complete ness obtained in o ur study with Montcgriffo's 
met hod described in IT'osi et all 120011) and the dashed lines are 
from Aloisi ct al. 1 2001]) (hereafter AL01). In the middle panels we 
show the photometric error distributions derived by AL01, whilst 
the distributions obtained with Montegriffo's method are shown in 
the bottom panels. We indicate as thick lines the cuts in the error 
distributions adopted to perform the simulations presented in this 
paper. 



derived by AL01 are plotted in Figures ^ The limit- 
ing magnitude is mpiiow — 24.5, and the color range 
is mFnow- m Fi60W = —0.8 - 3. To easily identify the 
main evolutionary p hases, we plotted the Padova tracks 
{Fagotto et all 11994ft for a metallicity of Z = 0.004 on 
the same scale (Fig. The dominant feature is the 
"red plume" extending from mpnow ~ 18 to 20.5. It 
has a small spread in color (0.2 mag) and a long range 
in magnitude. The region corresponds to the He burning 
phases of stars with masses approximately from 9 to 20 
M©. The more massive stars in our catalog have masses 
of 15 — 20 M©, corresponding to ages of 10 — 15 Myr. 
A second feature is defined by the bulk of objects with 
niFiiow ^21 — 24 and mFnow- m Fi60W > 0.9 where the 
tip of the RGB could be located. 

Unfortunately, the RGB feature is not recognizable in 
the observed CMD. The tip of the RGB should be lo- 
cated at niFiiow — 22.4 for an adopted distance mod- 
ulus DM = 26.71 and E(B - V) = 0.56. On the other 
hand, the DM uncertainty introduces an uncertainty of 
0.6 mag in the location of the RGB tip (between 21.8 
and 23.0 in mpnow)- I n addition, the RGB tip fea- 
ture on the LF could be blurred in the CMD because 
th e strong SF occurred at in termediate ages, as found 
bv lVallenari fc Bomansl l)1996|) and in this paper (see be- 
low). 

3. FROM THE EMPIRICAL TO THE SYNTHETIC CMDS 



The derivation of the SFH from the CMD is critically 
dependent on the completeness and analysis of the pho- 
tometric errors. This is particularly true for the oldest 
populations in a crowded field such as the central region 
of NGC 1569, where blending of unresolved stars may 
severely affect our data and create spurious features in 
the derived CMD. It is therefore crucial to quantify these 
effects as well as possible by following the recipes of the 
data reduction, and include them in the construction of 
the synthetic CMDs. 

3.1. Artificial star experiments 

We accounted only for the completeness of the F160W 
frame in our simulations because, as mentioned above, 
all the stars detected in F160W were previously found in 
F110W. However, we have performed the artificial star 
experiment (ASE) in both images, since the photometric 
errors in both bands are required. 

We initially adopted the ASE performed by AL01. 
They performed tests in which the fake stars have a PSF 
equal to that adopted in the data reduction. Then a frac- 
tion (~10%) of the number of detected stars was added 
at random positions on the image in each run and about 
10 runs were performed. The total number of artificial 
stars was then similar to that of the detected objects. 
They derived an effective limiting magnitude of the data 
of mpi60W — 23 and a 50% completeness level at mpi6ow 
~ 20. In the top panels of Figure [3 we plot their com- 
pleteness and photometric error distributions. The re- 
sults obtained for mpigow (mFiiow) are shown in the 
right (left) panel. In the top panels we plot the com- 
pleteness factors from AL01 as dashed lines, while in the 
middle panels we draw their estimated photometric error 
distributions. 

When we applied the ASE results from AL01 to the 
simulations, we found that they do not allow us to re- 
produce the faintest part of the CMD. For instance, if 
we simulate an old SF episode (starting 5 x 10 9 yr ago 
and still active) with the same number of stars of the 
observed CMD, we obtain a synthetic diagram that un- 
derproduces the faint (niFiiow ~ 23.5) tail of the LF by 
50% . We interpret this result as being due to an under- 
estimate of the true completeness in the AL01 ASE, and 
we therefore performed new ASE computations. We be- 
lieve that this is because ALOl's approach does not add a 
sufficient number of artificial stars at faint magnitudes, 
where crowding, blending and incompleteness are very 
severe. We therefore adopted a procedure de veloped by 
P. M ontegriffo at the Bologna Observatory l|Tosi et alJ 
I2001() . which uses ~ 10 5 artificial stars per band. We 
have put fake stars in the images according to a LF sim- 
ilar to the observed one. At brighter magnitudes this 
LF follows the observed one, but at low luminosities, a 
higher number of stars is required because of the high 
probability of losing objects because of crowding. All 
the fake stars have a PSF equal to that adopted in the 
actual data reduction. To avoid artificial crowding, we 
divided the frames in grids of cells of known width, and 
we randomly positioned only one fake star per cell at each 
run. We forced each star to have a distance from the cell 
edges large enough to guarantee that all its flux and back- 
ground measuring regions fell within the cell. In this way, 
we could control the minimum distance between adjacent 
stars. At each run, the position of the grid is randomly 
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changed and the artificial stars are then searched in the 
image by using the same procedure adopted for the real 
data and selected with the same criteria. 

An artificial star was considered lost when the differ- 
ence between the input and output magnitude was larger 
than 0.75. Such a difference implies that the artificial 
star overlaps with another star of the same luminosity or 
brighter. In the top panels of Figure|2 the completeness 
derived using Montegriffo's method is shown as a solid 
line. The photometric error distributions are plotted in 
the bottom panels as dots. The effective limiting magni- 
tude of the data is mpi60W — 24.5 and the sample is 50% 
complete at mpi6ow — 21.5. We found a posteriori that 
the completeness function obtained with Montegriffo's 
method allows us to reproduce the faint LF significantly 
better than with the previous ASE. 

The input-output magnitude difference of the artifi- 
cial stars provides an estimate of the photometric errors, 
which is safer than the error indicators of the reduction 
package {o-Daovhot hi our case). The outcomes are plot- 
ted in panels |2b and c for the ASE approaches by AL01 
and Montcgriffo, respectively. The synthetic CMDs ob- 
tained by adopting the photometric errors from the arti- 
ficial star tests, usually reproduce the color and magni- 
tude spread of the obser v ational CMDs very well [se e 
e.g. iGallart et all 120021). iSmecker-Hane et all (2002), 
iTosi et all l|2002fl . lAnnibali et all l|200.^ ]. In this case, 
however, by adopting the input-output magnitude differ- 
ence of Fig. [2t> as photometric error, we obtain an exces- 
sively large color spread in the sequences of the synthetic 
CMD. This effect is clearly evident along the red plume, 
where the observed CMD is tight in color. We therefore 
used the synthetic CMDs to heuristically infer the ap- 
propriate error distribution in each filter. First, we cut 
the wings of the photometric error distributions at the 
la level in each magnitude bin, and adopted these dis- 
tributions to derive a test CMD. Since in this test CMD 
the ler photometric error still has a color distribution 
larger than observed, we iterated the cutting procedure 
until we obtained photometric error distributions capa- 
ble of reproducing the spread in the main features of the 
CMD. The cuts performed on the photometric error dis- 
tributions eventually adopted are those shown within the 
solid lines in the bottom panels of Fig. 

3.2. Preliminary tests 

We have also checked whether different portions of the 
image have different completeness and photometric er- 
rors. To this end, we have divided the NIC2 images into 
three regions: the central one covering the most crowded 
part of the image including the SSCs, and the other two 
towards the edges in opposite directions. No relevant de- 
viation has been found either in the completeness factors 
or in the input-output magnitude difference of the three 
regions. 

Many authors have found evidence of spatially 
variable intrinsic extinction in NGC 1569; this ex- 
tinction increases toward the central region of the 
galaxy, where SSCs and H I I regions are located 
llGonzalez-Delgado et al.l Il997t iDeyost. Rov &: Drissenl 
Il997t iKobuhiickv fc Skillmanlll997HOriglia et al.ll2001fl . 
We have divided the image in subregions and compared 
the corresponding color functions in order to investigate 
whether strong differential reddening is present in our 



field. No significant differences in the color distributions 
have been detected from the image center to the edges. 
We are confident that differential reddening is not an 
issue in our 205 x 205 pc 2 sized region. 

Since the CMDs of all the subregions show similar pho- 
tometric errors, completenesses and color distributions, 
we have derived a global SFH over the whole 205 x 205 
pc 2 image. 

3.3. The synthetic CMD method 

The synthetic CMDs are constructed via Monte Carlo 
extractions of (mass, age) pairs, according to an as- 
sumed IMF, SF law, and initial and final epoch of the 
SF activity. Each synthetic star is placed in the theo- 
retical log{L / Lo),log(T e f f) plane by suitable interpola- 
tions on the adop ted evolutionary tra cks (in this case, 
the Padova sets bv lFagotto et al.lli994j) . Luminosity and 
effective temperature are transformed into the desired 
photometric system by interpolation within appropriate 
tables for photometric conve rsions [the transformation s 
to the VEGAMAG system bv lOriglia fc Leithererl |2000), 
in this case] , and the resulting absolute magnitude is con- 
verted to a provisional apparent magnitude by applying 
reddening and distance modulus. An incompleteness test 
on the F160W magnitudes is then performed, based on 
the results from Montegriffo's ASE. A photometric error 
is assigned to all the retained stars by using the "mod- 
ified distributions" discussed in Section 3.1. The value 
of the SFR is obtained when the number of objects in 
the synthetic CMD (or in regions corresponding to spe- 
cific evolutionary phases) is consistent with the observa- 
tional data. When the examined CMD region contains 
only a few objects, we perform several simulations in or- 
der to avoid problems due to small number statistics. In 
these cases, we adopted the average value of the resulting 
SFRs. 

4. SIMULATIONS 

In order to derive the SFH of a galaxy, we explore a pa- 
rameter space dependent on the distance modulus, red- 
dening, IMF, SF law, metallicity, and age of the episodes. 
We obtain a solution when we have a combination of pa- 
rameters that generate a synthetic CMD similar to the 
observed one. The exploration is also limited by the in- 
trinsic uncertainties of the theoretical models, like those 
introduced by the stellar evolution tracks and the pho- 
tometric conversions. 

First, we attempt to reproduce the observed CMD with 
a single episode and metallicity model over the whole 
range of ages covered by the empirical CMD. If a sin- 
gle episode does not reproduce the observed CMD, we 
assume multiple SF events. If a single metallicity does 
not allow us to simulate an acceptable CMD, we try to 
change the metallicity for each episode. Then we per- 
form other simulations in order to check the effect on 
the resulting SFH by assuming different IMFs, distance 
moduli and reddenings. 

This approach gives us an indication of the uncertainty 
on the results of our simulations. It allows us to iden- 
tify the range of the parameter values leading to a sat- 
isfactory agreement with the data, and to reject those 
clearly inconsistent with the observational evidence. Fol- 
lowing this line of thought, we will show the simulated 
episodes that fit the observed data but also some of the 
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Fig. 3. — Top panels SIM1: an old single-episode. The galaxy 
forms stars from 5 X 10 9 yr ago until today. The mean SFR is 0.28 
Mq yr -1 kpc - 2 . Middle panels SIM2: a young single-episode. In 
this case NGC 1569 forms stars from 5 X 10 s yr ago to 1 X 10 7 
yr ago with a mean SFR ~ 1.1 Mq yr -1 kpc -2 . Bottom panels 
SIM3: a single-episode that forms stars from 1.5 X 10 s yr ago to 
1 X 10 7 yr ago with a SFR of 1.7 Mq yr -1 kpc -2 . 



bad episodes, because they provide hints for a better un- 
derstanding of the SFH. We ran more than 1000 simula- 
tions to infer the SFH of the complicated case of NGC 
1569. 

We first present the results obtained by adopting 
DM = 26.71, E(B - V) = 0.56, Z = 0.004 and a Salpeter 
IMF (a =2.35) over the whole mass range from 0.1 to 120 
M Q . In paragraph 4.3 we discuss the effect of assuming 
different values for the IMF slope, distance and redden- 
ing parameters. The models all assume a constant SFR 
within the specific age ranges; we always report the SFR 
of each episode both in M© yr -1 kpc -2 and in M Q yr -1 
(the latter in brackets). 

4.1. Single episode SF 

We have run several cases with a single episode of 
SF in NGC 1569, all constrained to reproduce the total 
number of stars in the observed CMD. Here we describe 
only three representative models: the first one (hereafter 
SIM1, plotted in Fig.|3|) assumes a long SF activity, com- 
mencing 5 Gyr ago and still ongoing. The second model 
(hereafter SIM2, plotted in Fig. |3J| assumes a SF com- 
mencing only 500 Myr ago and terminating 10 Myr ago. 
The third simulation (hereafter SIM3, plotted in Fig- EJ 
generates stars from 150 to 10 Myr ago. 

SIM1 has a SFR of 0.28 M© yr" 1 kpc" 2 (0.01 
Mq yr -1 ). The simulated objects are concentrated in the 



lower, fainter part of the CMD (mpnow ~ 22.5), where 
the low-mass stars in the RGB phase are located. The 
observed CMD does not show a similar concentration 
of objects in the faintest part. Conversely, at brighter 
magnitudes (mpiiow < 22.5), the synthetic red plume 
shows a star deficiency with respect to the observed LF. 
These inconsistencies show that the bulk of the assumed 
SF activity is located at excessively early epochs and pro- 
vides too many old stars. At the same time the synthetic 
CMD is populated by blue massive stars (with M > 20 
M Q ). A dozen of these synthetic stars are in the main- 
sequence phase. However, we cannot clearly distinguish 
the main-sequence from the other evolved phases in the 
observed CMD because of the large photometric errors 
at faint magnitudes. We have also performed several 
simulations to overcome the small number statistics of 
the stars located at bright luminosity (mpnow < 18.5). 
These simulations suggest that too many red supergiant 
stars are produced without observational counterparts. 
We can conclude that, in addition to the excess of old 
stars in the synthetic CMD, we also have an excess of 
young stars. This is the reason why in SIM2 we have 
turned on the SF more recently and turned it off earlier 
(10 Myr ago). 

SIM2 has a SFR - 1.1 M yr" 1 kpc" 2 (~ 0.04 
Mq yr -1 ) and fails to account for the brighter stars of 
the red plume. Again, at faint magnitudes, SIM2 shows 
an excess of stars. In both simulations (SIM1 and SIM2), 
we have a feature in the region 19.5 < mpiiow < 21.5 
and 1.45 < mFnow- m Fi60W < 1-7, which corresponds to 
stars with masses 4 — 6 M© in the AGB phase and does 
not appear in the observed CMD. We then tried to sim- 
ulate the CMD with a single episode starting 1.5 x 10 8 
yr ago and terminating 10 Myr ago in order to avoid 
these stars. SIM3 has a SFR of 1.7 M Q yr" 1 kpc" 2 (0.07 
M© yr" 1 ). This simulation avoids the feature obtained 
in SIM1 and SIM2 but shows a shortage of stars at low 
luminosities (mpiiow > 23) and at the top of the red 
plume (mpnow < 20), and an excess of objects in the 
middle (20 < mpiiow < 23) of the luminosity function. 

At the end of these preliminary simulations, no solu- 
tion has been found that simultaneously fits both the 
bright and the faint parts of the CMD. Hence, a single 
episode cannot reproduce the observed CMD: to improve 
the fit to the data, we need to perform multiple-episode 
cases. 

4.2. Multiple episodes SF 

We use the age-box procedure to investigate multiple 
episodes. The basic idea is that we can perform our 
analysis by dividing the CMD in regions (boxes) popu- 
lated by stars with selected prop erties, such as age or 
evolutionary p hase iGreggidl2002l) . A similar procedure 
is adopted bv lGallart et alJ l)1996|) . The star counts in 
each box are directly linked to the mass that went into 
stars in the SF episodes that occurred at the correspond- 
ing epoc hs, i.e. the ave rage SFR over the specific time 
intervals l)Greggioll2002fl . 

In order to better understand the age-box procedure, 
we introduce Figure 0] where we plot the apparent mag- 
nitude versus age of red (mFiiow _ni Fi60W > 0.9) post 
main sequence stars. Vertical shading indicates to He- 
burning stars, and slanted shading indicates RGB (angle 
= +45°) and E-AGB (angle = -45) stars. Because of 
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Fig. 4. — Magnitude range covered by red post-main-sequence 
stars of a simple stellar population with an age between 10 Myr and 
10 Gyr selected in color mpnow-iiiFi60W > 0-9- Vertical shading 
refers to He-burning stars; slanted shading to RGB (angle = +45) 
and E-AGB (angle = -45) stars. This plot clearly shows that the 
temporal resolution is worse at older ages. 



the evolutionary lifetimes, the red stars at magnitudes 
brighter than 22.4 (where the tip of the RGB is located) 
are mostly core-He burning objects. The straight relation 
between their magnitude and age (Fig. 0J thick shaded 
strip) translates into a good temporal resolution of the 
SFH at ages younger than about 0.1 Gyr. This allows 
us to obtain precise results in the boxes corresponding 
to the top of the red plume. At magnitudes fainter than 
22.4, each LF bin collects (in principle) stars with ages 
up to the Hubble time. The stellar counts in this region 
allow us to derive only integrated information. 

In practice, we first define the magnitude ranges of 
the box by identifying the features in the LF 4 that sug- 
gest a possible discontinuity in the star formation. The 
color ranges are chosen in such a way to account for the 
spread because of the photometric errors. These ranges 
constrain the box that defines our partial model, which 
is constrained to reproduce the observed stellar counts 
in the box with a SF activity within a specific age range. 
The limits for the age range are read from Fig. 01 in 
correspondence to the magnitude limits of the box. Be- 
sides the objects in the box, the simulation also populates 
the diagram in regions different from the box described 
above; we take these "additional" objects into account 
when we perform the simulations in the other parts of 
the CMD. We have run hundreds of cases to infer the 
best combination of parameters. However, for the sake 
of clarity, here we simply describe the line of reasoning 

4 The comparison between model and data is performed selecting 
the red LF. However, in this paper we report only on the total 
luminosity function, since our CMD is mostly composed by a red 
population and the red LF does not significantly differ from the 
total one. 



and show the most relevant results from all these tests. 
Th is method has been app li ed an d extensively d i scusse d 
bv ISchulte-Ladbeck et ail llgOOOT): ICrone et alJ (12001 : 
iSchulte-Ladbeck et all <|2002D : IAnnibali et all <|2003[) . 

4.2.1. One box simulations 

Following the above reasoning, we define the box (here- 
after BOX1) in the magnitude range 18.5 < mpnow 
< 20.5 (at mpiiow = 20.5 a dip in the LF is notice- 
able) and colors 1.0 < mFnow _ni Fi60W < 1-4 (to take 
into account the photometric errors). Figure 01 shows 
that in BOX1 we are sampling core-He burning stars in 
the age range (1 — 4) x 10 7 yr; these objects are the 
youngest, most massive stars of the red plume. This box 
also includes a population of lower-mass stars in the AGB 
phase; because of their short lifetimes in this magnitude 
range, these stars provide a minor contribution to the 
overall population. 

After many simulations, we find the best parameter 
combination (SIM4) that reproduces the number counts 
in BOX1 (Fig. EJ). SIM4 is an episode commencing 
3.7 x 10 7 yr ago and finishing 1 .3 x 10 7 yr ago with a mean 
SFR of 3.1 - 3.3 M yr" 1 kpc" 2 (0.12 - 0.14 M yr" 1 ). 
We performed other simulations by slightly changing the 
boundary ages of the episodes for an estimate of the un- 
certainty affecting the determination of the age intervals. 

Other possible solutions have a starting epoch in the 
range (3.5 — 3.8) x 10 7 yr ago, and an ending epoch in the 
range (1.3 — 1.4) x 10 7 yr ago. In all the explored cases, 
the simulations for this young episode share a problem 
that is also visible in SIM4: an overproduction of stars 
with mFiiow- m Fi60W < 0.5 (hereafter blue bulk). This 
structure corresponds to the blue edge of the blue loop 
phase of 8 — 12 M Q . We cannot remove this feature with- 
out underproducing the objects in the red plume. The 
number of objects generated in this recent SF episode is 
20% of the total. 

We have then tried to simulate the remaining objects 
with a single episode, but we have not succeeded in re- 
producing the LF. If we construct single episodes start- 
ing at epochs earlier than 3 x 10 8 yr ago and terminating 
4x 10 7 yr ago, we obtain too many faint objects (mpnow 
> 22.5) and too few stars in the range mpnow 20.5 — 22.5 
in the synthetic diagram. On the other hand, an episode 
that begins to form stars ^3x 10 s yr ago, does not pro- 
duce a sufficient number of objects at faint magnitudes. 

The lack of success in reproducing the LF with only 
two episodes led us to consider other scenarios. 

4.2.2. Two box simulations 

The second box (hereafter BOX2) is defined in the 
magnitude range 20.5 < mpnow < 22.5 and color 0.9 < 
niFnow- m Fi60W < 1-3 and targets the intermediate ages. 
The bright limit is fixed by BOX1, while the faint limit is 
chosen to be slightly brighter than the RGB tip in order 
to exclude possible blended RGB stars from the counts. 
From Fig. 0J one can see that stars falling in BOX2 have 
ages in the range 1.2 x 10 8 yr to 4 x 10 7 yr, if they are 
core-He burners. The color limits of BOX2 were indeed 
chosen with the aim of isolating core-He burning stars, 
while taking into account the color spread due to photo- 
metric errors. 

The best simulation for BOX2 (hereafter SIM5) ac- 
counts for ~ 50% of the objects and is plotted in the 
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Fig. 5. — The best synthetic CMDs and luminosity functions. 
In the top panel the youngest episode (SIM4) that reproduces the 
objects in BOX1 is plotted. The galaxy experienced strong SF 
beginning 3.7 X 10 7 yr ago and terminating 1.3 X 10 7 yr ago with 
a high rate of ~ 3.2 Mq yr - 1 kpc -2 . In the second panel from 
the top we plotted the intermediate episode (SIM5) constrained by 
BOX2. This episode occurred in the interval 1.5 X 10 8 — 4 X 10 7 yr 
ago with a SFR ~ 1.2 Mq yr" 1 kpc -2 . The oldest episode that 
reproduces the remnant objects (SIM6), is reported in the third 
panel from the top. Age interval: 6 X 10 s — 3 X 10 8 yr ago. 
SFR = 1.5 M yr -1 kpc -2 . The final composed synthetic CMD 
is drawn in the bottom panel. 



second panels from the top in Fig. [3] SIM5 has a SFR = 
1.2 M yr -1 kpc -2 (0.04 M Q yr -1 ). This episode started 
forming stars 1.5 x 10 8 yr ago and terminated 4.0 x 10 7 
yr ago. The starting epoch is older than the 1.2 x 10 s 
yr mentioned above. When computing the simulations, 
we vary the starting and ending epoch around the values 
derived from Fig. 0] in order to fit the observed LF with 
the minimum number of individual SF episodes. It turns 
out that starting at 1.5 x 10 s yr allows us to fit the stellar 
counts in the range 22.5 — 23.5 without adopting a new 
box. Notice that at such magnitudes the completeness 
factors fall below 50%. 

If we assume that a quiescent phase (or one with low 
SFR) occurred (3.7 — 4) x 10 7 yr ago, we can also repro- 
duce the shape of the LF in the magnitude range mpiiow 



20 — 20.5. This kind of short inter-burst phases has also 
been found by G98. 

4.2.3. The remaining objects 

The simulations computed to reproduce the stellar 
counts in BOX1 and BOX2 account for about 70% of 
the total number of observed objects. The remaining 
30% must be ascribed to SF that occurred at epochs 
older than 1.5 x 10 s yr. From Fig. 0] we see that the 
CMD in the magnitude range 22.5 — 24.5 can be popu- 
lated with a mixture of relatively young He-burning stars 
[with ages of (1.5 — 3) x 10 s yr], early- AGB stars with 
intermediate ages (up to 1 xlO 9 yr), and RGB stars with 
ages up to the Hubble time. Therefore, we tested various 
starting epochs for the oldest episode of SF sampled by 
our data. 

First, we attempted to account for the remaining stars 
by considering He-burning objects with ages between 
3 x 10 8 yr and 1.5 x 10 s yr. The simulation fits the LF 
fainter than mpnow < 22.5 reasonably well, but yields 
too many stars in the region 19.5 < mpnow < 21 at col- 
ors 1.45 < mFnow _m Fi60W < 1-7. The overproduction 
amounts to a factor of 4 — 5. These are bright E-AGB 
stars with masses between 4 and 6 M . As discussed in 
Section 4.1, these synthetic objects were also produced 
in SIM2, which assumed active SF at intermediate ages. 
If we use the number counts in the box 19.5 < mpnow 
< 21; 1.45 < mFiiow- m Fi60W < 1-7 to constrain the SF 
activity at ages between 1.5 x 10 8 and 3 x 10 8 yr ago, 
we obtain a SFR level of 0.14 M yr -1 kpc" 2 (0.006 
Mq yr -1 ), which is 1 order of magnitude lower than 
the rate of NGC 1569 in the young and intermediate 
episodes. We then tested SF episodes setting in at var- 
ious epochs, but all ending 3 x 10 8 yr ago, to avoid the 
overproduction of these bright AGB stars. We found that 
if the old SF episode started a few Gyr ago, the upper 
RGB stands out very clearly on the synthetic CMD, with 
a prominent horizontal feature corresponding to the tip. 
This is not seen in the data. For a starting epoch of 
2 x 10 9 yr ago, the signature of the RGB tip is no longer 
apparent, and we obtain consistency with both the CMD 
and the LF. 

We obtain acceptable models for any starting epoch in 
the range 2 - 0.6 Gyr, with corresponding variations of 
the SFR. For example, a simulation started 1.5 x 10 9 yr 
ago yields a SFR of 0.4 M yr -1 kpc" 2 (0.009 M Q yr -1 ); 
a simulation started 6 x 10 8 yr ago yields SFR = 1.5 
M Q yr -1 kpc -2 (0.06 M yr -1 ). Our best fit is actu- 
ally the latter, since the former case tends to underpro- 
duce the faint AGB stars, corresponding to 3 Mq objects 
(20.5 >mFiiow > 21.5) . The best simulation for the re- 
maining objects is shown in Fig. [S] as SIM6, while the 
best composite CMD is plotted in Fig. [S] with the label 
Tot. 

4.3. Searching for other parameter combinations 

We have performed many other simulations in order to 
search for alternative solutions; in particular, we consid- 
ered cases with different IMFs, with slope a in the inter- 
val 1.35 — 3.35. Fig.[fj]plots the SFR as a function of the 
IMF exponent, for the youngest (circles) and intermedi- 
ate (triangles) episodes. Fig. also shows that the SFR 
has a minimum at a m i n ~ 1.5 — 2. This behavior is due to 
the non-linear relation between the slope of the IMF and 
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Fig. 6. — The trend of the SFR as a function of the IMF slope 
(a). Filled circles refer to the youngest episode (BOX1, 8 Mq < 
M < 15 Mq) and filled triangles to the intermediate one (BOX2, 
5 Mq < M < 8 Mq). Dashed lines show the theoretical relation 
between the IMF slope and the SFR of the corresponding box. The 
differences are due to stochastic fluctuations. 



the number of stars per unit mass of the parent stellar 
population. By adopting a steep IMF (e.g. a ~ 3), the 
low-mass stars are relevant in number and weight with 
respect to the total population, while the contribution of 
high-mass stars to the total mass is negligible. For a flat 
IMF (e.g. a ~ 1.3), low-mass stars are still numerous, 
but their weight decreases. On the other hand, the num- 
ber of massive objects become slightly higher but their 
contribution to the total mass increases more than the 
total number. 

We can broadly reproduce the shape of the LF in 
BOX1 with an IMF flatter than Salpeter's (down to 
a = 1.35). For simplicity, in this paper we have adopted 
a single-slope IMF, extending Salpeter's IMF to the mass 
range 0.1 - 120 Mq. IGould. Bah call fc Flvnnl l)1997|) sug- 
gest a slope natter tha n Salp eter's (a — —0.56) for M 
< 0.6 M Q . iChabrierl l|2003l) has also found that the 
IMF in the Galaxy depends weakly on the environment 
and is well described by a power-law for M > 1 Mq 
and a log-normal form below this mass, except pos- 
sibly for early SF condit ions. For the Galactic disk, 
iKroupa fc Weidnerl l)2003ft found an IMF with a slope 
a = 1.3 in the range 0.08—0.5 Mq and a slope that is con- 
sistent with Salpeter's above 0.5 Mq. If we extrapolate 
the SFR at M < 0.5 Mq by adopting Kroupa's rather 
than Salpeter's IMF, the values become 30% smaller than 
with the default Salpeter IMF. 

We have performed quantitative comparison (x 2 ) be- 
tween our best model and alternative solutions without 
the gap in the age range 37—40 Myr, and found that the 
model with the gap is ~ 2000 times more probable than 
the fit without the gap. 

Another possible SFH reproducing the brightest por- 



tion of the red plume consists of many short episodes 
in the age interval (1 — 4) x 10 7 yr ago. This possi- 
bility is also indicated by the "granularity" of the LF. 
In this case we obtain the best fit with three episodes 
separated by short periods (few million years) of qui- 
escence or very low SFR. SIM4a started 2.1 x 10 7 yr 
ago and terminated 1.3 x 10 7 yr ago with a SFR of 
~ 4.2 M Q yr" 1 kpc~ 2 (~ 0.18 M Q yr" 1 ); SIM4b started 
3.0 x 10 7 yr ago and terminated 2.2 x 10 7 yr ago with a 
SFR of ~ 4.5 M Q yr" 1 kpc" 2 (~ 0.19 M Q yr" 1 ); SIM4c 
started 3.7 x 10 7 yr ago and terminated 3.3 x 10 7 yr ago 
with a SFR of ~ 3.2 M Q yr" 1 kpc" 2 (~ 0.13 M yr" 1 ). 
This SFH reproduces remarkably well the brightest por- 
tion of the total LF. These three episodes have SFRs 
equal or higher than the single episode (see SIM4) but 
their total astrated mass is similar to that. We prefer not 
to pursue this fine tuning (although we cannot exclude 
a more complex scenario in the recent star formation) 
because the shape of the red plume and its LF may be 
affected by residual star-like systems mistaken as indi- 
vidual stars in spite of all our checks (see AL01) . 

There are no clear evolutionary features in the ob- 
served CMD (Fig. ^) that allow a precise estimate of 
the distance modulus of the galaxy. Therefore, we have 
performed several simulations in order to examine the 
effects of changing the distance on the derived SFH in 
NGC 1569. To this end, we have adopted the lowest 
and the highest values of the DM found in the literature 
(DM = 26.1 — 27.3). In general, we obtain a different 
fit to the data by changing the DM. For sake of clar- 
ity, we only discuss the case of the larger distance, since 
this could be important in order to understand whether 
the RGB tip can be hidden. In this case, we derive a 
younger episode with higher SFR with respect to our 
best model. This is because a given apparent magnitude 
samples intrinsically brighter, and hence more massive, 
stars when DM is larger. As a consequence, a larger 
number of massive and young stars needs to be repro- 
duced by the simulations. Also, when we modify the 
distance of the galaxy, we also change the position of the 
evolutionary features of the synthetic CMD. This effect 
is important because it could substantially modify the 
SFH. By adopting DM = 27.3, the best fit is obtained 
with a young episode that has formed objects in the age 
interval 2.7 x 10 7 — 8 x 10 6 yr ago, and the intermediate 
SF occurred 1 x 10 8 — 3.2 x 10 7 yr ago. The correspond- 
ing SFRs are about 2 times higher than those found for 
DM = 26.71 in SIM4 and SIM5. A SF with a rate of 0.8 
M Q yr -1 kpc~ 2 (0.04 M yr -1 ) that occurred in the age 
interval 1.3 x 10 9 — 2 x 10 8 yr ago is consistent with the 
data. We can exclude a strong SF occurring at ages older 
than 1 — 2 Gyr because the RGB tip would be clearly 
present in the CMD despite the higher distance. We want 
to stress again that, for what concerns the older SF, the 
results must be qualified by considering the uncertain- 
ties at faint magnitudes because of the incompleteness 
and photometric errors. 

The reddening uncertainty suggested by the spread of 
the literature values shifts the synthetic CMD by 0.05 in 
color (mFnow- m Fi60w) and 0.1 in magnitude (mpiiow); 
its effect on the age of the episodes is smaller than those 
due to the uncertainties on the distance modulus and 
tracks. 

Finally, with the synthetic CMD method, we can pro- 
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vide hints on the metallicity of the SF episode. Our 
capability to select the metallicity of an episode depends 
on the photometric errors relative to its main evolution- 
ary phase. In principle, we are able to reject all the 
metallicities that create unobserved features in the syn- 
thetic CMD. The colors of the red plume are best re- 
produced with Z=0.004; the Z=0.0004 tracks yield red 
supergiants that are too blue compared to the obser- 
vations. The Z=0.008 tracks instead produce very red 
objects. At fainter magnitudes, where the photometric 
errors are larger, the presence of a low Z population can- 
not be ruled out, while the Z=0.008 tracks produce stars 
that are too red with respect to the data. 

4.4. Observational and theoretical uncertainties 

The low galactic latitude of NGC 1569 (b ~ 11°. 2) led 
us to check the number of Milky Way stars falling into 
ou r field of view. AL01 ra n the Galaxy model described 
bv lCasertano et aD l|1990(l and found that the foreground 
contamination in our (small) analyzed region is negligi- 
ble. Major problems in our analysis are instead crowding, 
the shallow limiting magnitude and the high reddening. 
Crowding is severe in the central region because the high 
recent SF activity has generated many bright, young ob- 
jects that hide fainter, older ones. In addition, the shal- 
low limiting magnitude does not allow us to reliably de- 
tect features that could help us to better determine the 
characteristics of the older episode. From this point of 
view, despite the high reddening, deeper, higher resolu- 
tion HST images in optical bands could have possibly 
provided a better data set to study the earlier epochs. 
This h as been clearly demonstrated by lAnnibali et al.l 
( 2003) in NGC 1705 (a BCD twice as distant as NGC 
1569). In that case, WFPC2 photometry has allowed to 
reach 1 — 1.5 mag below the RGB tip in regions where 
NIC2 barely allows the identification of the tip. 

An additional uncertainty on the analysis of the central 
region of NGC 1569 comes from t he presence of many 
small clusters {Hunter et all 120001: lAnders et all 12004(1 
and star-like systems. Most of them appear as bona fide 
stellar objects in our photometry, with shape, sharpness, 
and x 2 parameters typical of single stars. These unre- 
solved systems could be located in the CMD at the top 
of the red plume. By keeping or removing all possible 
cluster candidates from our empirical CMD, we estimate 
that the star-like systems could change the SFR by 15 — 
20% and shift the end epoch of the youngest episode to 
an age a few Myr older. 

As discussed above and shown in the top panels of 
Fig. EI what we called blue bulk in Sect. 4.2 could affect 
the derivation of the SFR of the youngest episode. This 
spurious feature could be due either to uncertainties in 
the bolometric corrections in the infrared filters or to 
uncertainties in the stellar tracks at the blue edge of the 
blue loops for the masses 8 — 12 M© . In order to evaluate 
the influence of the blue bulk, we re-simulatcd the top of 
the red plume by locating a box in the region of the blue 
bulk (20.5 < mpnow < 21.5 and 0.9 < mFnow- m Fi60W 
< 1.3). By using the same age interval of the youngest 
episode, the simulation generates a partial model, plot- 
ted in Fig.0 which provides a SFR of 1.6 M Q yr -1 kpc -2 
(0.06 M© yr" 1 ). This simulation is able to generate only 
50% of the total objects in the red plume of the observed 
CMD and is clearly inconsistent with the data. We also 
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Fig. 7. — Simulation attempting to reproduce the number of 
objects in the box of the blue bulk. The number of objects at the 
top of the red plume is underestimated: this is an indication that 
the SFR found in SIM4 is an upper limit. 



performed many simulations by using Padova tracks of 
different metallicity, but we found similar problems. In 
other words, either we reproduce the number of stars 
in the blue bulk region, but severely underproduce the 
red plume objects, or we reproduce the red plume and 
overestimate the blue bulk. We consider the latter op- 
tion more viable, since the timescales of stellar evolution 
models at the blue edge of the blue loops are more un- 
certain than those o f the red plume phases. Notice that 
IGreggio et alJ ljl998f) have derived the SFH in the B and 
V bands adopting both the Geneva and Padova tracks 
in order to evaluate the uncertainty related to different 
stellar evolutionary models. They found an excess of 
synthetic stars at the blue edge of the blue loop with the 
Geneva tracks but not with the Padova tracks. This sug- 
gests that an additional uncertainty comes into play in 
the blue bulk] for instance, the completeness factor could 
be overestimated because its color dependence has not 
been considered in the simulations. 

5. DISCUSSION AND CONCLUSIONS 

The best agreement between the simulations and the 
data is obtained by adopting the SFH plotted in Figure|Hl 
(solid line). Keeping in mind that these results are af- 
fected by many uncertainties (both in data and models), 
we can derive the following overall scheme for the SFH 
of NGC 1569 (Recall that the results described in this 
paper refer to the central region of the galaxy). 

- NGC 1569 has experienced a complex SFH, com- 
posed of at least three strong episodes in the last 1 
(possibly 2) Gyr. 

- The quiescent phase that occurred (1.5 — 3) x 10 8 
yr ago is inferred from the small number of 3 — 6 
M© stars in the AGB phase. 

- We need an episode older than 3 x 10 s yr to re- 
produce the fainter part of the CMD. This old SF 
may start as early as 1 — 2 Gyr ago. The large 
uncertainties (severe incompleteness, large photo- 
metric errors) that characterize the fainter part of 
the CMD do not allow us to infer reliable details 
on the SF of the oldest episode. However, we find 
hints for a significantly lower SF activity (or even a 
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Fig. 8.— The derived SFH of NGC 1569. The SFH over the 
Hubble time is reported in the main plot linearly in age. The 
solid lines refer to the episodes of Figure [F] and the dashed line 
to the second scenario described in the conclusions. The recent 
star formation history resulting from the simulations is shown in 
the small logarithmic plot. We give only an upper limit for the 
youngest episode. These uncertainties in the SFR are due to the 
blue bulk feature. Dotted lines indicate other possible solutions 
described in the text. 



quiescent phase) at ages older than 1 — 2 Gyr (for 
a distance of 2.2 Mpc). 

A Salpeter IMF reproduces the brighter part of the 
observed CMD and the LF reasonably well. A flat- 
ter IMF is also able to reproduce the brighter por- 
tion of the red plume. 

By assuming a Salpeter IMF over the whole mass 
range, the most recent SF occurs at a rate of 3.1 — 
3.3 Mq yr" 1 kpc~ 2 (0.12 - 0.14 M© yr" 1 ), the 
intermediate one at a rate of 1 — 1.2 M© yr -1 kpc -2 
(0.03 — 0.05 M Q yr^ 1 ), and the older one at a rate 
between 0.25 and 1.5 M Q yr" 1 kpc" 2 (0.01 - 0.06 
M© yr -1 ), depending on the assumed duration. 

The uncertainties of the distance modulus affect 
the estimate of the ages and SFRs of the episodes. 
The best fit obtained with a larger distance of 
DM = 27.3 is 2.7 x 10 7 -8 x 10 6 yr ago for the young 
and 1 x 10 8 — 3.2 x 10 7 yr ago for the intermedi- 
ate episode. The corresponding SFRs are about 2 
times higher than those found for DM = 26.71. 
A SF with a rate of 0.8 M Q yr" 1 kpc" 2 (0.04 
M© yr -1 ) occurred in the age interval 1.3 x 10 9 — 
2 x 10 8 yr ago. Also for the larger distance, we can 
exclude a strong SF that occurred at ages older 
than 1 — 2 Gyr because the RGB-Tip would be 
clearly present in the CMD. 

We cannot exclude more than one SF episodes oc- 
curring in the age interval (1 — 4) x 10 7 yr . A pos- 



sible scenario could consist of three episodes with 
a SFR of 4.2, 4.5 and 3.2 M yr" 1 kpc~ 2 (0.18, 
0.19 and 0.13 M yr -1 ) separated by short phases 
of quiescence or low activity. The uncertainties are 
due to the possible presence of clusters at the top 
of the red plume and small number statistics. 

By comparing our results and those from the litera- 
ture for the SFH of NGC 1569, we find overall agree- 
ment, although with interesting differences. The end of 
the youngest episode (1.3 x 10 7 yr ago) is significantly 
older than that obtained bv lVallenari fc Bomansl (1996) 
and G98, who place it around (5 — 10) x 10 6 yr ago, as 
found from optical HST data. These differences can be 
explained in terms of a different sensitivity of the differ- 
ent bands (optical vs. NIR) to hot and cold evolution- 
ary phases. Here we are studying young massive stars 
with NIR data. We are then sensitive to the post-main- 
sequence phases, where these stars have very short life- 
times. Optical data, more appropriate for main-sequence 
objects, allow for higher reliability on the youngest stars. 
On the other hand, the NIR analysis allows us to dis- 
tinguish the two recent episodes, separated by a short 
quiescent period, better than with the optical data. In 
addition, we find evidence for a gap, or a low SF from 
3 x 10 s yr ago to 1.5 x 10 8 yr ago , similar to that suggested 
bv lVallenari fc Bomansl l)1996f) . 

The SFR derived by G98 using the PC camera (which 
covers an area of 374 x 374 pc 2 for a distance of 2.2 Mpc), 
is 4 M© yr -1 kpc~ 2 in the age interval 1.5 x 10 8 — 5 x 10 6 
yr ago. From Fig. [5J we can see that our rate for the 
youngest episode (SIM4) is slightly lower than that de- 
rived by G98. For the intermediate episode, we have 
derived a SFR 3 — 4 times lower than G98. We per- 
formed further simulations by adopting the ingredients 
suggested by the literature cited above in order to bet- 
ter understand these differences. We cannot reproduce 
the top of the red plume if we simulate the totality of 
objects in the CMD with a continuous and constant ac- 
tivity starting (1 — 1.5) x 10 8 yr ago and ending 5 x 10 6 
yr ago. Furthermore, this simulation underestimates the 
fainter part of the LF. Another simulation (SIM7, plot- 
ted in Fig. EJ) assumes the SFR of G98 in the same age 
interval. The resulting synthetic diagram has twice as 
many stars than our observed CMD. The total luminos- 
ity function is broadly fitted in the brighter part of the 
red plume but shows an overabundance of objects in the 
range mpiiow = 19-5 — 23.5, whereas the fainter part is 
barely populated. From the comparison of the observed 
LF with that from SIM7, we can learn that: 

- If we terminate the SF 5 x 10 6 yr ago, we produce 
many high-mass stars without observed counter- 
parts; 

- In the BOX1 area, SIM7 produces a number of 
objects ~ 1.5 times higher than SIM4. Neverthe- 
less, the SFRs of SIM7 and SIM4 are quite similar 
within the uncertainties since the length of the two 
simulations is slightly different. 

- The simulated objects of SIM7 in the BOX2 region 
are about 3 — 4 times higher than SIM5 (this value 
is similar to the ratio of the corresponding SFRs). 
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Fig. 9. — Simulation of an episode with duration and SFR as 
found by G98. These parameters do not allow us to reproduce the 
observed CMD and LF. 



This discrepancy could arise from the uncertain- 
ties in the photometric error distribution and com- 
pleteness factors used in this work and in G98. In 
addition, the differences in the SFH could result 
from the differences in the areas sampled by the 
two cameras. 

By using archival HST data. lAnders et all l)2004f) find a 
high-level of star cluster formation for ages younger than 
2.5 x 10 7 yr and a lower level of star cluster formation at 
older epochs, with a secondary peak about 100 Myr ago. 
The cluster age distribution shows a gap in the age inter- 
val (1.6-4.0) x 10 8 yr. These results match our SFH that 
refers to field stellar population. However. I Anders et al.1 
l)2004[) find a younger age for the end of the last episode , 
consistent with the results bv lVallenari fc Bomansll)1996f > 
and G98. This again suggests that the optical data are 
more appropriate for studying the most recent epochs. 

In order to estimate whether the expanding superbub- 
bles can be plaus ibly driven by the obse rved population 
of massive stars, iHeckman et al.l l)1995l) found that the 
observations can be reproduced reasonably well by as- 
suming an episode with a constant SF over the ages 
1.3 x 10 7 — 3.2 x 10 7 yr ago. The ages of the episode 
are rather similar to our youngest burst. 

By using popu lation synthesis models, 
iBelczvnski et alJ 1)20041) have constructed synthetic 
X-ray binary populations for direct comparison with the 
X-ray LF and found two principal stellar populations. 
One is old, metal-poor, with continuous SF from 1.5 
Gyr ago and the other is recent and metal-rich. Their 
results are consistent with ours for the oldest episode, 
but invoke a recent burst younger than that inferred 
here. Again, we attribute our earlier termination to the 
selection effect which makes it difficult to detect blue 
young stars in the NIR. We remind the reader that the 
most recent termination epoch for the SF inferred in 
the optical by G98 only refers to the regions of NGC 
1569 where single stars were resolved. There are zones 
in this galaxy occupied by giant HII region complexes 
where even with the HST/WFPC2 we were unable to 
resolve individual stars. The existence of these HII 
regions clearly suggests that the SF activity must have 
continued until even more recent times. 

All three episodes have SFRs higher by 2-3 or- 
ders of magnitude than that of the solar neighbor- 



hood [e.g. iTimmes. Wooslev &; Weaver! ()1995D ] and than 
those inferred with similar methods in other late-type 
dwarfs [both irregular s and BCDs, s ee the revi ews by 
ISchulte-Ladbeck et all pOOl and bv lTosil ^OO^i }. The 
only dwarf analyzed so far with the synthetic CMD 
method which has a SFR comparable to that of NGC 
1569 is NGC 1705 l|Annibali et al.l 1200^) . These two 
dwarfs, in spite of being classified differently (NGC 1705 
is a BCD), have quite similar properties, all somehow re- 
lated to the strong SF activity. They both contain SSCs, 
show evidence of galactic wi nds, and have SFRs con- 
sistent with those required bv lBabul fc Ferguson! 
to let dwarfs account for the excess of blue galaxies at 
intermediate redshift. Neither shows evidence of long 
quiescent phases in the last couple of Gyrs. 

At the rates derived in this paper, the central region of 
NGC 1569 has formed ~ 2.7 x 10 7 M Q of stars in t he last 
1 Gyr i.e. about 8% of the total mass estimated bv llsraell 
( 1988) and about 14% of the mass that is not associated 
with neutral hydrogen (M dif = 2 x 10 8 M Q ). The field of 
view of NIC2 corresponds to an area of 0.04 kpc 2 , which 
cover s only ~ 2.5% of the optical size of NGC 1569 [1.76 
kpc 2 l|Tsraelll1 988ft ]. 

To improve our knowledge of the SF history we have 
compared the integrated optical and infrared fluxes of 
the whole field and of the resolved stars in the NICMOS 
area, from the images described by G98 an snd AL01. 
The integrated magnitudes in the HST F439W, F555W, 
F110W and F160W filters are 13.7, 12.9, 10.9 and 10, re- 
spectively, for the field; 5 they are 15.2, 14.5, 12.1, and 11 
for the resolved stars. The dereddened mF439w — rriFi60W 
color of our investigated region is then ~ 1.6, noticeably 
bluer t han the corre sponding color for a solar Kurucz 
model l)Kuruczlll992l. (mF439w — m Fi60w)o — 2.2). The 
apparent integrated B magnitude of our NlCMOS field 
is mr^NiCMOS — 13.6, whil st the total B magnitude o f 
NGC 1569 is m B = 11.86 l|De Vaucouleurs et alJll991^ . 
Hence, our NICMOS frame samples ~ 20% of the total 
B light from the galaxy, if the average reddening of the 
NICMOS region is representative of that of the whole 
galaxy. Notice that our sampled area is only 2.5% of the 
total, therefore the SF per unit area in NGC 1569 is very 
centrally concentrated. 

We can explore some possible scenarios in order to 
place constraints on the past history of NGC 1569. Given 
the to tal and HI masses of NGC 1569 estimated bv llsraell 
( 1988) and the astrated mass in the last 1 Gyr extracted 
from our CMD procedure and assuming a Salpeter's 
IMF over the whole mass range, we are left with only 
1.73 x 10 s M Q for the whole galaxy for the sum of the 
older stars and the dark matter. Assuming that these 
old stars were produced at a constant rate over the in- 
terval 1.3 x 10 10 — 1 x 10 9 yr ago, the corresponding 
SFR is 1.44 x 10~ 2 M yr" 1 . Since this SFR refers 
to the whole galaxy, the normalized value is 8.2 x 10~ 3 
M Q yr -1 kpc -2 . This value is comparable to what has 
been derived via simulations for dwarf galaxies in the Lo- 
cal Group (see G98 and references therein) . For the NIC 
area, the corresponding astrated mass over the age inter- 
val is Mnic — 4.3 x 10 6 Mq. If we simulate an episode 
that generates the M^ic mass over the age interval re- 

5 These values do not include the flux from the SSCs and from 
the HII regions, which were masked on the frames. 
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ported above, we obtain a synthetic CMD that contains 
5 - 10% of the total objects of the observed CMD. These 
synthetic stars are on the RGB. Hence, in this scenario 
we cannot rule out that a constant but low SF might 
have occurred over the Hubble time. 

However, it would be very unusual for a galaxy 
of this type to have a dark/visible mass ratio much 
lower than 10. Moreover, chemi cal evolution models 
(R omano. Tosi fe Matteucci 1 I2005|) show that without 
the usual dark matter ratio, NGC 1569 would experi- 
ence too early strong galactic winds and it would be 
impossible to reconcile the model predictions with the 
observational data. 

In the second scenario (Fig. 8 dashed line) we there- 
fore assume that the stellar component is only 10% of the 
dynamical mass; then our astrated masses are probably 
inconsistent with a Salpeter's IMF extrapolated down to 
0.1 M Q . However if we adopt iKrouoa fc Weidnel ( 2003) 
IMF, the central region has generated only ~ 60% of the 
stellar mass of the galaxy over the last 1 Gyr. The re- 
sulting mean SFR over the interval 1.3 x 10 10 - 1 x 10 9 
yr ago is 6.7 x 10~ 4 M Q yr _1 kpc~ 2 . Alternatively, con- 
sidering an onset of the SF 7 Gyr ago (i.e., the galaxy 
started to form stars at redshift 1 by adopting the "con- 
cordance cosmology" [Hq = 70, fl\ = 0.7, flu = 0.3]), 
the mean SFR is about 2 times higher than the value 
above. All these quantities lie in the range of SFRs typ- 
ical for nearby dwarf galaxies. These scenarios do not 
allow us to rule out a past SF with a low rate (or similar 
to the typical SFR of the other dwarf galaxies), but we 
can safely conclude that the last 1 — 2 Gyr have been 
peculiar. 
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What could be the cause of this peculiar SF in the last 
Gyr? Interestingly, there is an HI cloud with mass 7 x 10 6 
M Q located at a projected distance of 5 kpc from NGC 
1569. It is connected to the galaxy by a bridge similar in 
mass to the HI cloud ijStil & Israel 19 93)- This cloud ma y 
have triggered the star formation (Stil & Israel 1998). 
Moreover, according to Muehle (2003), who studied the 
galaxy with high-resolution HI data, the structure of the 
galaxy halo is the remnant of an intergalactic HI cloud 
in the process of merging with the galaxy. If true, this 
may provide both the cause and the fuel for the strong 
and concentrated SF activity of NGC 1569. 

To better understand the history of this intriguing ob- 
ject, we should try to resolve its stars older than the few 
Gyr reached by our current data, in both the inner and 
outer regions. Although the NIR bands are important to 
explore old stellar population, NICMOS is not an opti- 
mal instrument. The detection of the tip of the RGB (if 
it exists) will be achievable only with the HST Advanced 
Camera for Surveys, which has both the necessary spatial 
coverage and resolution. 
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